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Using phase-equivalent supersymmetric partner potentials, a general result from the inverse prob- 
lem in quantum scattering theory is illustrated, i.e., that bound-state properties cannot be extracted 
from the phase shifts of a single partial wave, as a matter of principle. In particular, recent R-matrix 
analyses of the 12 C + a system, extracting the asymptotic normalization constant of the 2 + sub- 
threshold state, C12, from the I = 2 elastic-scattering phase shifts and bound-state energy, are 
shown to be unreliable. In contrast, this important constant in nuclear astrophysics can be deduced 
from the simultaneous analysis of the I = 0, 2, 4, 6 partial waves in a simplified potential model. A 
new supersymmetric inversion potential and existing models give C12 = 144.5 ± 8.5 x 10 3 fm" 

PACS numbers: 03.65.Nk,25.55.Ci,26.20.+f 
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I. INTRODUCTION 

Can bound-state properties be deduced from scatter- 
ing data? The answer to this question showed interesting 
evolutions throughout the history of quantum scattering 
theory (see Foreword by R. G. Newton in Ref. start- 
ing with Heisenberg's conjecture in the early forties that 
the scattering (S) matrix contains all physical informa- 
tion, including that regarding bound states. This con- 
jecture was later seriously weakened in the framework 
of the potential model: in 1949, Bargmann 0, con- 
structed potentials that share the same S matrix for par- 
tial wave I = (so-called phase-equivalent potentials be- 
cause they share the same scattering phase shifts) but 
have different bound-state properties [different energies, 
or identical energies but different wave functions]. On 
the other hand, that same year, Lcvinson (sec for in- 
stance Ref. pj) proved that, for a given partial wave £, 
the number of bound states, Ni, can be deduced from 
the phase shifts. These results were then explained by 
the inverse-problem theory which shows that for a 
given partial wave the potential is uniquely determined 
by (i) the scattering phase shifts at all energies, (ii) the 
bound-state energies and (iii) one additional real param- 
eter for each bound-state wave function [e.g., the asymp- 
totic normalization constant (ANC), see definition JIJ]. 
This only applies to local central potentials displaying 
no singularity. Levinson's theorem, however, was gener- 
alized by Swan in 1963 0, who proved that N# cannot be 
deduced from the S matrix for potentials displaying an 
r~ 2 repulsive singularity at the origin. Since then, the 
construction of phase-equivalent potentials for a given 
partial wave has known important progress 0, 0, 0, IS 13 
thanks to the use of the algebraic formalism of supersym- 
metric quantum mechanics [lJJ, 0, ^3 used below (see 
Fig.0. 

As far as only one partial wave is concerned, the answer 
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to the above question is thus rather clear-cut: neither 
the number of bound states, nor their energies nor their 
ANCs can be deduced from the S matrix alone. More- 
over, these bound-state properties are independent of one 
another, which means that knowing, for instance, both 
the S matrix and a bound-state energy does not constrain 
the value of the ANC as long as no other physical infor- 
mation is available on this particular bound state (e.g., 
its radius or its lifetime) or on the considered system in 
general (e.g., the range of the interaction). This impos- 
sibility proof being established in the framework of the 
potential model, which is a particular case of more gen- 
eral models (e.g., microscopic or phenomenological mod- 
els), it holds for these more general models too. However, 
several works 0, Q, 0, 0, 0] seem to ignore this ar- 
gument and aim at extracting the ANC of a bound-state 
with known energy from the corresponding partial-wave 
S matrix. In Refs. for instance, the reaction- (R- 

)matrix formalism |l8|,ll9j * s used to extract the ANC of 
an 16 bound state from 12 C + a elastic-scattering data. 
In the R-matrix phenomenological model, it is rather nat- 
ural to assume that bound states have a measurable im- 
pact on the partial-wave S matrix, like elastic-scattering 
resonances, because both bound and resonant states are 
characterized by an energy and a reduced width. How- 
ever, in scattering theory, bound and resonant states have 
very different natures, and the above impossibility ar- 
gument seriously questions the ANC extraction of Refs. 
|l6l Il7| . In Sec. [ft] I show how that misuse of the R- 
matrix formalism actually explains the inconsistency be- 
tween the ANC values obtained in these references. Then 
in Sec. IIII1 I explicitly construct simple potential mod- 
els which illustrate the independence of the bound-state 
ANC from the phase shifts for a single partial wave. 

Now, considering several partial waves simultaneously 
leads to a less clear situation. The inverse-problem 
theory implies that a regular interaction potential is 
uniquely determined by phase shifts of a given partial 
wave with angular momentum high enough so that the 
partial wave contains no bound state. Using inversion 
techniques E|, a unique potential can thus in principle 
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be deduced from such phase shifts; this potential can in 
turn be used for lower partial waves, for which it may 
sustain bound states. In that case, bound state prop- 
erties (of low partial waves) are totally determined by 
scattering phase shifts (of a high partial wave). How- 
ever, this program presents three difficulties. First, it 
assumes that the potential is local and energy indepen- 
dent, which, in general, is not a realistic approximation 
when the interacting particles have an internal structure 
(see the discussion in Sees. El an d [3 for the particular 
case of 12 C + a). Second, it requires that the interac- 
tion be independent of the partial wave, which is not 
always the case: scattering phase shifts of several partial 
waves are not necessarily compatible with one unique 
potential. Third, it is based on the knowledge of the 
phase shifts of a high-angular-momentum partial wave 
at all energies, which are generally not available: at high 
energy, the complexity of multichannel effects increases, 
which makes the data difficult to analyze; and at low en- 
ergy, the phase shifts of high-angular-momentum partial 
waves are very small because of the centrifugal barrier. 
Hence, considering several partial waves at the same time 
may help deducing bound-state properties from scatter- 
ing data, but the practical applicability of the method 
could strongly depend on the considered system. In Sec. 
IIVI I show that for the particular case of 12 C + a encour- 
aging results are obtained by considering several partial 
waves at the same time, in the framework of a simplified 
model in which 12 C and a are considered as rigid nu- 
clei in their + ground state, interacting through a local 
energy- and angular-momentum-independent potential. 
These results are then compared with more sophisticated 
models available in the literature. Conclusions and per- 
spectives for the 12 C + a system are given in Sec. Ivl 



II. R-MATRIX ANALYSES 

The 12 C + a system is particularly important in nu- 
clear astrophysics since the 12 C (a, 7) 16 capture is 
a key reaction in the helium-burning phase of red giant 
stars: by competing with the triple a reaction (leading to 
12 C), it determines the ratio of carbon and oxygen which 
results from stellar nucleosynthesis. Experimentally, the 
Coulomb repulsion between 12 C and a makes the direct 
measurement of the capture reaction impossible at the 
very low energies of astrophysical interest (around 300 
keV in the center-of-mass frame) . Hence, available exper- 
imental results at higher energies (-E c-m , > 1 MeV) have 
to be extrapolated to these low energies. This extrap- 
olation, generally performed with the R-matrix model 
(the only model to date able to precisely fit all the avail- 
able data), is itself problematic because of the presence 
of two 16 bound states just below the 12 C -I- a thresh- 
old. Whereas the influence of the Ec.m. = —45.1 ± 0.1 
keV state is well understood [20I l2lj . the influence of 
the -E c . m . = —244.9 ±0.6 keV state is still the major 
source of uncertainty on the reaction rate [l6l ITtI 211. 



The unknown key quantity is its ANC (the energies of 
both states are well known), which is related to its a 
reduced width in the R-matrix formalism (I do not use 
reduced widths here since they depend on the R-matrix 
radius, which makes comparison between different mod- 
els difficult). 

The —245 keV subthreshold state has a positive parity 
and a J = 2 angular momentum (2 + state); microscopic 
cluster models (see Ref. [2^ and references therein) sug- 
gest that its structure is dominated by a 12 C cluster and 
an a cluster, both in their fundamental + state, with a 
relative motion of angular momentum 1 — 2. It is thus 
on the phase shifts of this partial wave that one could 
naively expect to see the influence of the 2 + state, and 
it is from these phase shifts that one could try to extract 
its ANC, C12 (where the second index refers to I — 2 
and the first index means that this state is the lowest 
one for this partial wave). This constant is defined by 
the asymptotic behavior of the radial part, u\i(r), of the 
relative wave function between the 12 C and a nuclei for 
this particular bound state: 

«i 2 (r) ~ C 12 W_ 6 . 615 . 5/2 (0.3751 r), (1) 

r — >oo 

where r is the radial coordinate (in fm) and W is the 
Whittaker function [23} (a decreasing exponential func- 
tion in the absence of Coulomb interaction). 

In Ref. ^(| , the 1 — 2 phase shifts from Ref. [24| are an- 
alyzed by the R-matrix model, leading to a minimum in 
the x 2 (see Fig. 1 of Ref. 0) for C 12 = 402 x 10 3 fm- 1 / 2 . 
However, this extraction is very delicate since the influ- 
ence of the subthreshold state has to be disentangled from 
the background phase shifts in the energy region where 
experimental data are available; this influence is smaller 
(see Fig. 2 of Ref. 0|) than the main components of this 
background, i.e., the hard-sphere and high-energy-pole 
phase shifts. As a consequence, this minimum strongly 
depends on the range of data fitted (see Fig. 1 of Ref. 
|l6j). Moreover, I have checked that the minimum is 
highly sensitive to the high-energy background-pole en- 
ergy, which is fixed at 10 MeV in Ref. [l6j. Without 
this restriction, very different values of C\i are obtained; 
hence, the error on C12, which is not estimated in Ref. 
(l6j . is actually very large. 

In Ref. 01 1 the 12 C + a elastic scattering has been re- 
measured with very high precision. An R-matrix analysis 
with free high-energy-pole energies leads to a x 2 mini- 
mum (see Fig. 2 (a) of Ref. [HI) for C 12 = 154 ± 18 x 10 3 
fm" 1 / 2 . The error bar on C12 is estimated by following 
the guideline [U HI 

+ 9Xmh>, (2) 

where Xmin ^ s the minimal x 2 arL d v is the number 
of degrees of freedom. In the present c ase, one has 
X 2 ^ = 18941 (the x 2 of Fig. 2 in Ref. Hfl have to 
be divided by the 32 reference angles), for 11392 data 
points and 65 parameters (32 R-matrix parameters and 
33 experimental parameters), which leads to v = 11327. 
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Following the theory of x 2 minimization [25|. one finds 
that the probability of having x 2 > 18941 for v = 11327 
is of the order of 10 -391 , which probably indicates that 
some systematic errors have been underestimated in Ref. 
[l7| . Consequently, Eq. cannot be used to define error 
bars on parameters since this equation only takes into ac- 
count statistical errors and requires that the fit be good 
in the sense of x 2 theory. A less rigorous way of checking 
the validity of a fit is to use the "chi-by-eye" approach 
|25|. In the case of Ref. 0], the fit looks good (see 
Fig. 1 of this reference), which is related to the fact that 
Xmin/ 17 = 1-672 is not that far from 1. However, the fits 
obtained with C u = 136 or 172 x 10 3 fin -1 / 2 , which cor- 
respond to x 2 l v — 1-674, look probably equally good. I 
thus consider that the error bar on Cyi obtained in Ref. 
[l7j is underestimated. 

In summary, the values of C12 obtained in both Refs. 
[l6| and 0] are most probably unreliable, for different 
reasons. This explains why, despite the fact that the 
I = 2 phase shifts of Refs. and [13, [gjj are com- 
patible with one another, their R-matrix analyses give 
totally incompatible results for C\i. One might think 
that this inconsistency is due to the fact that in Ref. 0] 
all relevant partial waves are taken into account whereas 
in Ref. [16| only £ — 2 is considered. However, since the 
R-matrix formalism is based on a partial-wave decompo- 
sition of the cross section, C\i is actually independent of 
the £ ^ 2 parameters (in Ref. 0], the only influence of 
the £ ^ 2 phase shifts on £ = 2 comes through the choice 
of a common R-matrix radius for all partial waves; how- 
ever, this does not seem to strongly constrain the fit of the 
£ = 2 phase shifts, as discussed in the next paragraph). I 
argue that this incompatibility rather illustrates the im- 
possibility, well-known in quantum scattering theory, of 
deducing a bound-state ANC when only the bound-state 
energy and the elastic-scattering phase shifts of the cor- 
responding partial wave are known. 

Let us however notice that, in the case of the R-matrix 
phenomenological model, this impossibility is probably 
not as strong as in general scattering theory. As shown 
in the next Section (see Fig.[2J), the range of a 12 C + a in- 
teraction that reproduces both the £ = 2 phase shifts and 
the bound-state energy increases with the value of C\2- 
But the R-matrix formalism is only valid when the nu- 
clear interaction is negligible above the R-matrix radius. 
In the R-matrix theory, where the number of R-matrix 
poles is infinite [n|, this does not cause any problem: 
the results are independent of the R-matrix radius, which 
can thus be chosen as large as necessary (in other words, 
it would be possible to calculate an R-matrix for any 
of the potentials plotted in Fig. [2] but the R-matrix ra- 
dius should be larger than the potential range). In the 
R-matrix phenomenological model, where the number of 
poles is limited to the number of bound states and reso- 
nances in the energy interval considered, plus one possi- 
ble high-energy pole to help simulating the background, 
a constraint appears on the R-matrix radius: in Ref. [27| . 
a value of 5.5-6 fm is found, based on the £ = 1 partial 
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FIG. 1: Experimental I = 2 background phase shifts de- 
duced from the 12 C + a elastic-scattering data of Ref. [rHl26| 
(points), compared with the theoretical phase shifts of the 
phase-equivalent potentials of Fig. [5] (solid line, identical for 
all potentials). 

wave; in Ref. 0] [see Fig. 2 (b)] 5.5 fm leads to the best 
fit (I have checked that this constraint is mainly due to 
the £ = partial wave, for which the background is max- 
imal); in Ref. ^(|, the radius is fixed at 6.5 fm. Hence, 
very large values of C12 can probably be excluded by an 
R-matrix fit; very small values, on the other hand, can- 
not be excluded. I could actually find very good fits of 
the £ — 2 phase shifts with C12 = for R-matrix radii in 
the range 4.5-6.5 fm. This confirms the doubts expressed 
above on the analyses of Refs. [rIITtT ] and it shows that, 
in this particular case, the quality of the fit does not 
strongly depend on the precise value of the R-matrix ra- 
dius. 



III. INVERSION OF ONE PARTIAL WAVE 

Up to now, the discussion is based on general argu- 
ments from scattering theory. Let me now construct 
simplified two-body potential models that illustrate these 
arguments in the particular case of 12 C + a, using the 
supersymmetric-quantum-mechanics [Tol ITlT Wl\ inver- 
sion technique of Refs. |28|, |29j. Starting from a phase- 
shift analysis 1261 o f the purely elastic data (E c , m , < 4.966 
MeV) of Ref. [T3, I have removed, using the R-matrix 
formalism, resonances that have a more complicated 
structure than 12 C + a in their + ground states (such 
resonances cannot be explained by the present simplified 
model). For £ — 2, this means removing the narrow reso- 
nance at 2.683 MeV and the wide resonance at 4.358 MeV 
[30| . which leads to the phase shifts plotted in Fig. ^ (the 
spread of the data above 4.2 MeV is due to the difficulty 
of performing a reliable phase-shift analysis in the vicin- 
ity of a resonance). By solving the so-called "singular 
inverse problem" defined in Ref. [28j], I have constructed 
a potential that reproduces these phase shifts with very 
high precision. This potential has an r~ 2 singularity at 
the origin and satisfies Swan's theorem |j; to simplify 
further developments, its nuclear part has been approxi- 



4 




1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 



,,20 




1 1 








/ \ r 




2000 















2 4 6 8 10 12 
r (fm) 



FIG. 2: Effective 12 C + a interaction potentials for I = 2. 
The solid-line potential corresponds to Eq. © while the 
dashed-line potentials are its phase-equivalent supersymmet- 
ric partners with one bound state at —245 keV, for C12 = 20, 
200 and 2 000 x 10 3 fm" 
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FIG. 3: Radial part of the 2+ subthreshold state 12 C + 
a relative wave functions. The dashed lines correspond to 
the potentials of Fig. |5] the solid line to the deep inversion 
potential and the dash-dotted line to the potential of Ref. 



mated by the analytical expression 

V(r) = Fexp(-r 2 /i? 2 )/r 2 , (3) 

where V = 43.40 MeV fm 2 and R = 5.091 fm. The cor- 
responding £ = 2 effective potential (centrifugal barrier 
plus Coulomb and nuclear interactions) is represented in 
Fig. |21 (solid line) and its phase shifts are shown in Fig. 
n Surprisingly, a potential of the form is also able 
to reproduce the £ — background phase shifts with 
satisfactory precision, which suggests an alternative way 
of parameterizing the 12 C -I- a background in future R- 
matrix analyses (in the spirit of the "hybrid" R-matrix 
model of Refs. |3lLl32.l33l |V The fit of the £ = 2 scattering 
phase shifts is very good, despite the fact that potential 
@ has, by construction, no bound state (it is purely re- 
pulsive); this is a first proof that, from the point of view 
of scattering theory, the £ = 2 background phase shifts 
are totally disconnected from the 2 + subthreshold bound 
state. 

Now, since this state is known to exist and to have 
a structure such that it could be well simulated by the 
present simplified model, it is physically meaningful to 
perform a further step in the inversion of the £ — 2 
phase shifts, i.e., to add a bound state to the above 
potential without modifying its phase shifts |2^. This 
phase-equivalent bound-state addition may be carried 
out through a further use of supersymmetric quantum 
mechanics 0; it introduces two arbitrary parameters: 
the bound-state energy and its ANC. In the present case, 
the energy is known experimentally but not the ANC; 
hence, there is a one-parameter family of £ = 2 effective 
potentials that share identical phase shifts (Fig. ^) and 
have a bound-state at —245 keV. Such potentials are rep- 
resented by dashed lines for three arbitrary values of C\i 
in Fig. |2 they constitute generalizations (to £ ^ and 
in the presence of a Coulomb interaction) of Bargmann's 
potentials (see Fig. 1 in Ref. 0). The corresponding 
bound-state wave functions are plotted in Fig. [3] (dashed 
lines), where their very different asymptotic behaviors 



are clearly seen. This proves that, from the £ — 2 phase 
shifts and bound-state energy only, there is no way of 
telling what is the best value for C12, in contrast with 
R-matrix attempts described above (see however the dis- 
cussion of the R-matrix phenomenological model in the 
previous Section). 



IV. INVERSION OF SEVERAL PARTIAL 
WAVES 

The general discussion from the introduction suggests 
to look at other partial waves simultaneously. The 2 + 
subthreshold state is actually a member of a rotational 
band, together with the + bound state at —1 113±1 keV 
and the 4 + and 6 + resonances at 3 194 ± 3 keV (width: 
26 ± 3 keV) and 9 113 ± 7 keV (width: 420 ± 20 keV) 
respectively [3(j. This is indicated by cluster-model cal- 
culations (see Ref. [2^ and references therein) , which give 
a fair description of these four states in a single-channel 
calculation where both the 12 C and a nuclei are in their 
+ ground state. This suggests to ask the following ques- 
tion: could all these states be reproduced by a unique 
^-independent potential? Of course, in principle, the in- 
ternal structure of the nuclei could play a role and a sat- 
isfactory potential could be non local and energy depen- 
dent. It is however worth testing the simplest hypothesis, 
i.e., that of an energy- and ^-independent local potential, 
because if such a potential were sufficient, it would be 
strongly constrained by the data of the four partial waves 
considered simultaneously, which could help solving the 
C12 indefiniteness. 

I could not find such a potential within the above one- 
parameter potential family. This is not very surpris- 
ing since general arguments from the microscopic cluster 
model show [35| that ^-independent interaction poten- 
tials between nuclei should have low-energy non-physical 
bound states that simulate the Pauli principle between 
the constituting nucleons [because of these so-called Pauli 
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forbidden states (PFSs), such potentials are deep, in con- 
trast with shallow potentials that only sustain physical 
bound states]. For the 1 = 2 partial wave of 12 C + a, 
the expected number of PFSs is two |3^| . 

In the context of inversion, such bound states com- 
plicate things a lot since each of them introduces two 
additional arbitrary parameters (neither their energy nor 
their ANC is known experimentally). Following Ref. [29| . 
I introduce these states phenomenologically by regular- 
izing the above inversion potential with the expression 
V Teg +ar +br for r < r reg , where V reg and r reg are fitting 
parameters (a and b are calculated to make the poten- 
tial and its first derivative continuous in r rcg ). Since this 
regularization slightly alters the phase shifts, parameters 
V and R in Eq. ((3J) have to be adapted to maintain the 
fit of Fig. n] Altogether, there are thus five parameters 
(V, R, Ci2, V rcg and r rcg ) which can be varied to fit all 
the £ = 0, 2, 4, 6 data at the same time. The best fit 
I have obtained is for V = 63.81 MeV fm 2 , R = 4.925 
fm, V reg = -107.9241 MeV, r rog = 4.1 fm. By construc- 
tion, this inversion potential perfectly reproduces the 2 + 
state energy. The wave function is represented in Fig. 
(solid line); its two nodes are due to its orthogonality 
to the PFSs at —46 and —18 MeV, which confirms that 
the potential is deep. The corresponding value for Cu is 
137 x 10 3 fm" 1 / 2 . 

I consider this prediction as reliable because all the 
even phase shifts in the elastic region are also rather 
well fitted (for example, the 4 + resonance is at 3 196 
keV with width 22 keV). However, this potential under- 
estimates the + bound-state energy (—1696 keV), as 
well as the 6 + resonance energy (8 985 keV) and width 
(259 keV). These inaccuracies are probably due to the 
extreme simplicity of the present model (energy- and l- 
independent local potential) and it is worth comparing it 
with more sophisticated models of the literature, which 
better take into account many-body effects. The two 
non-local ^-dependent potentials of Ref. [23 , which also 
fit the physical properties of the members of the rota- 
tional band, lead to similar wave functions for the 2 + 
subthreshold state (2 nodes). The Gaussian potential 
provides C12 = 144 x 10 3 fm -1 / 2 while the Woods-Saxon 
potential provides 153 x 10 3 fm -1 / 2 . In Ref. 0, a sim- 
ple local potential is given, which precisely reproduces 
the energies of the members of the rotational band, as 
well as the resonance widths, provided a slight £ depen- 
dence is allowed. Although the number of PFSs is three 
for this potential, which adds a node in the subthrcshold- 
bound-state wave function (see dash-dotted line in Fig. 
|3J|, it provides a value of C\ 2 very similar to the inversion 
potential: 136 x 10 3 fm" 1 / 2 . 



V. CONCLUSION 

I have illustrated on 12 C + a two general results from 
the inverse problem in quantum scattering theory: (i) 
for a given partial wave, the scattering phase shifts, the 



number of bound states, their energies and their ANCs 
are independent of one another; (ii) when several partial 
waves can be described by the same theoretical model, 
these quantities may be related to one another. This is 
reasonably the case for 12 C + a, for which I have deduced 
the ANC value of the 2 + subthreshold state from the si- 
multaneous analysis of the £ = 0,2,4,6 partial waves, 
with the help of a local energy- and angular-momentum- 
independent potential constructed by supersymmetric in- 
version. 

Since this potential lies on rather restrictive hypotheses 
and is not able to exactly fit all the data, I have compared 
it with results from more sophisticated potential models 
available in the literature. All theoretical predictions lie 
in a rather limited range, C12 = 144.5 ±8.5 x 10 3 fm -1 / 2 , 
which confirms somehow their overall validity. The error 
bar can thus be considered as a sensible estimate of the 
theoretical uncertainty on C12 in a 12 C + a potential 
model. The fact that both local (present work and Ref. 
|34) 1 and non-local [33| potentials give consistent results 
is an indication that non-local effects are not very im- 
portant for this particular state. This confirms that its 
structure is dominated by rigid 12 C and a clusters with 
£ = 2 relative angular momentum. A similar smallness of 
non-local effects has been found in Ref. [36] for the 16 O 
+ a system. 

All the potential models discussed above neglect mul- 
tichannel effects, which means that the predicted value 
of C12 might be overestimated. This would be consis- 
tent with the experimental value obtained from transfer 
reactions H2, C 12 = 114 ± 10 x 10 3 fm" 1 / 2 , and with 
the theoretical estimate from the microscopic multicon- 
figuration cluster model of Ref. (H, C 12 = 134 x 10 3 
f m -!/ 2 Whereas compatible with the unreliable value 
from Ref. [l7j, the above theoretical estimate excludes 
that from Ref. 0] , as well as the cascade-transition value 
C12 = 228±p x 10 3 fm" 1 / 2 from Ref. A new mea- 
surement of the cascade transitions is planned at TRI- 
UMF, which hopefully should clarify the situation. In 
the meanwhile, R-matrix fits should reflect the uncer- 
tainty on C12. 

Let me end with a general discussion of possible theo- 
retical models for 12 C + a and similar systems relevant 
to nuclear astrophysics. Ideally, ab initio calculations 
should be performed; however, despite the success of the 
microscopic cluster model, such calculations are still too 
primitive to provide precise predictions of all available 
data. Hence, phenomenological models still have an im- 
portant role to play. Among them, the potential model is 
probably the most physical and intuitive one. However, 
simulating the full complexi ty o f the many-bodynuclear 
system requires a non-local [33l ] multichannel [38j treat- 
ment. This considerably complicates the model and a 
good quality fit of all available data with the potential 
model is not available to date. In contrast, the R-matrix 
model allows such a fit because it remains rather simple 
even in the presence of many channels. The price to pay 
for that is a large number of parameters and a very phe- 
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nomenological description, in particular of background 
terms. An interesting approach which combines the ad- 
vantages of both the potential and R-matrix models is 
that of "hybrid" models [3ll 133 , l33l | . The potentials con- 
structed in the present work renew the interest for these 
models. In particular, purely repulsive potentials of the 
form could be used as a new description of the back- 
ground, a possibility which will be explored in a future 
work. 
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